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Abstract. This paper studies the asymptotic behavior of eigenvalues of random abelian 
G-circulant matrices, that is, matrices whose structure is related to a finite abelian group 
G in a way that naturally generalizes the relationship between circulant matrices and 
cyclic groups. It is shown that, under mild conditions, when the size of the group G 
goes to infinity, the spectral measures of such random matrices approach a deterministic 
limit. Depending on some aspects of the structure of the groups, whether the matrices are 
constrained to be Hermitian, and a few details of the distributions of the matrix entries, 
the limit measure is either a (complex or real) Gaussian distribution or a mixture of two 
Gaussian distributions. 

1. Introduction 

Given a finite group G and a function / : G — > C, the matrix M = [fi a b~ 1 )] ab&G 
is called a G-circulant matrix by Diaconis [7J [8]. This generalizes the classical notion of 
circulant matrices, which arise as the special case in which G is a finite cyclic group. The 
action of such a matrix M on the vector space {g : G — > C} is as a convolution operator: 
for g : G — > C and a € G, 

(1.1) (Mg)(a) = Y^f{ab~ 1 )g{b) =: (f * g)(a). 

b&G 

This paper considers the asymptotic behavior of the spectra of random G-circulant ma- 
trices, or equivalently random convolution operators on G, when G is a large abelian group. 
(For the rest of this paper, G will always stand for a finite abelian group.) Such random ma- 
trices will be generated by picking the values f(a) independently, with or without imposing 
a constraint /(a" 1 ) = f(a) which is equivalent to insisting that the matrix M is Hermit- 
ian. This generalizes the study of random circulant matrices, whose theory has already 
been developed in [31 HI [U [151 E] among many other papers, with applications discussed in 
\10\ [TS] . The richer structure of arbitrary abelian groups relative to cyclic groups leads to 
the appearance of some interesting phenomena which do not occur for circulant matrices, 
or the more familiar setting of random matrices with independent entries. 

The prototypical situation (exemplified in Corollaries I3.2| 13.4) and I3.6[ and Theorems 
14.21 and 14.41 below) is that when the size of G grows the empirical spectral distribution of a 
(properly normalized) random G-circulant matrix M approaches a Gaussian distribution. 
When M is constrained to be Hermitian the limit will be a real Gaussian distribution; 
without such a constraint it will be a complex Gaussian distribution. These situations 
may be thought of as analogous to the semicircle law for Hermitian random matrices and 
circular law for non-Hermitian random matrices with independent entries, respectively. This 
behavior, which has already been observed for random circulant matrices in [31 [15], occurs 
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in particular if only a negligible fraction of the elements of G are of order 2, and also if 
every nonidentity element of G is of order 2. On the other hand, if neither of these is the 
case then more complicated limiting distributions occur which are mixtures of two Gaussian 
distributions (as in Theorems 14.11 and 14.31 below). 

Another perspective on these results, which is crucial in the proofs, is that they describe 
the distribution of values of random Fourier series on G. The supremum of such a random 
Fourier series is already a thoroughly studied quantity |11] 113] . In particular, results of 
Marcus and Pisier [13] include as special cases estimates of the spectral norms of random 
G-circulant matrices, as pointed out in Proposition 12.41 below. 

Section [2] below briefly reviews the facts about Fourier analysis on finite abelian groups 
which are used here and points out their immediate consequences for G-circulant matrices; 
some notation and conventions used in the remainder of the paper are established there. 
Section [3] investigates the spectra of some random G-circulant matrices whose entries are 
Gaussian random variables. The invariance properties of Gaussian random variables allow 
an easy detailed study to be undertaken which illuminates the general situation, in particular 
the role of the number of elements of order 2. Finally, Section [4] determines the asymptotic 
behavior of the spectrum for general entries with finite variances. 

The cases of G-circulant matrices with heavy-tailed entries, and of random G-circulant 
matrices when G is a nonabelian finite group, will be investigated in future work. 
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2. Some Fourier analysis and notation 

For a finite abelian group G, we denote by G the family of group homomorphisms x '■ 
G — > T, where T is the multiplicative group {z G C | \z\ = 1}. The elements of G are 
called characters of G; G is a group under the operation of pointwise multiplication. The 
multiplicative inverse of a character x is its pointwise complex conjugate x- From the 
homomorphism property it follows that for a G G and x € G, x{ a ~ l ) = x( a )- 

We denote by £ 2 (G) the space of functions / : G — > C equipped with the inner product 

a&G 

and 1 2 {G) is defined analogously. The Fourier transform of / G £ 2 (G) is the function 
/ G £ 2 (G) given by 

/(x) = (/,x) = E/( a ^( fl )- 

aeG 

This includes as special cases both the classical discrete Fourier transform (when G is 
cyclic) and the Walsh-Hadamard transform (when G is a product of cyclic groups of order 
2). The following lemma summarizes the most important fundamental facts about the 
Fourier transform for our purposes. 

Lemma 2.1. Let G be a finite abelian group with \G\ elements. 
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(1) The functions |— | \ G G} form an orthonormal basis of £ 2 (G). 

(2) The map f i— >• ^J^ f ^ s a linear isometry of £ 2 (G) onto £ 2 (G). 

(3) If f,g G i 2 (G), then for each \ G G, f * g(x) = f(x)d(x) (where the convolution 
f * g is defined in f 1 1 . 1 j) . 

Proof. (1) See Theorem 6 on j!6l p. 19]. 

(2) This follows easily from Proposition 7 on [161 P- 20] (which is a consequence of part 
©)• 

(3) This follows directly from the definitions by a straightforward computation. □ 

Observe that contained in Lemma l2.1[ fTj) is the fact that \G\ = \G\. 

We will need two additional facts about characters of finite abelian groups which are not 
as easily located in standard references. 

Lemma 2.2. The number of elements a £ G such that a 2 = 1 is equal to the number of 
characters x £ G such that x = X- 

Proof. For a G G, define 5 a ■ G — > C by 5 a (b) = 5 a ,b, where the latter is the Kronecker delta 
function, and observe that {5 a \ a G G} is an orthonormal basis of l 2 (G). Then 5 a (x) = xi a ) 
for each x G G. By Lemma l2.1t [2]). the number of a G G such that a 2 = 1 is equal to 

aeG 1 1 adG 1 1 a&G x< zQ 

which by Lemma 12 . 1 f )Tj) is equal to the number of x £ G such that x = X- ^ 

Lemma 12.21 says that G and G have equal numbers of elements of order 2. A much 
stronger fact is also true: G and G are isomorphic groups. However, this isomorphism is 
noncanonical, depends on the classification of finite abelian groups, and in any case is not 
useful here. 

Lemma 2.3. Let H be a subgroup of a finite abelian group H . Then each character on H 
extends to a character on G in precisely \G\ / \ H\ distinct ways. 

Proof. It is easy to check that restriction to H defines a homomorphism G — > H. Since 
each coset of this homomorphism's kernel has the same size, it suffices to prove that that it 
is surjective, or equivalently that each character on H extends to a character on G at all. 
For a proof of this fact see, e.g., [H p. 134]. □ 

From (|l.ip and Lemma I2.1I |3|) it follows that the Fourier transform diagonalizes G- 
circulant matrices. In particular, if M = [f(ab~ 1 )] a ^eG f° r / G ^ 2 {G), then the eigenvalues 
of M are precisely the values {/(x) I X £ C} of the Fourier transform of /, and the char- 
acters of G are eigenvectors of M. (For generalizations of these facts for nonabelian G, see 
[TIE]-) Observe that every G-circulant matrix is normal, but that M is Hermitian if and 
only if /(a -1 ) = /(a) for each a G G. 
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Given a family of random variables {Y a | a 6 G}, define the random function / € £ 2 (G) 
by /(o) = (We are avoiding using X to name random variables because of its 



typographical similarity to X-) The corresponding G-circulant matrix is the random matrix 
M = [^ab-^a i, & q- Its eigenvalues, indexed by x G G, are given by 

(2-1) A x = /( X ) = -^=^yaX(a), 

V l G l a6 G 

and the empirical spectral distribution of M is 

where <5 2 here denotes the point mass at z € C. 

The Fourier transform / is a random trigonometric polynomial on G, of the kind studied 
extensively by Marcus and Pisier |13j . From (|2.1|) it follows in particular that ||M|| = ||/|| , 
where the former norm is the spectral norm of M. The following result is thus a special 
case of [13} Theorem 1.4], which also applies to infinite compact abelian groups. 

Proposition 2.4. Suppose that {Y a \ a £ G} are independent (except possibly for a con- 
straint Y a -i = Y a for each a 6 G) and mean with finite second moments. Then 



( minE \Y a \] < E H M H < C.lmaxE\Y a 
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where c, C > are constants, independent of G and the distributions of the Y a . 

The rest of this paper deals mainly with infinite sequences of finite abelian groups G^ , 
always assumed to satisfy \G^ | — > oo. For each n a family of random variables {Yg | g £ 
G( n )| will be used to construct a random G^^-circulant matrix 



M (n) 



1 V (n) 

2 aft" 1 



with empirical spectral measure // ( - n ). As mentioned earlier, an important role will be played 
by the quantity 

(n) = |{a£GW|fl 2 = l}| = \{ X £G\x = x}\ 
P2 |G( n )| |G(")| 

The standard real Gaussian measure is denoted 7k, and the standard complex Gaussian 
distribution, normalized such that E|Z| =1 when Z is a standard complex Gaussian 
random variable, is denoted 7c- For a G [0, 1], j a denotes the Gaussian measure on C = M 2 
with covariance I f 1 ^" i so that in particular 70 = 7c and 71 = tr. 

The integral of a function / with respect to a measure ^ will be denoted by v(f). 

3. Gaussian matrix entries 

The following is an immediate consequence of Lemma l2.1Hl2j) and the rotation- invariance 
of the standard Gaussian distribution. The special case of this result for classical circulant 
matrices (that is, when G is a cyclic group) was observed in [15] . 



THE SPECTRA OF RANDOM ABELIAN G-CIRCULANT MATRICES 



5 



Proposition 3.1. Let G be a finite abelian group and let {Y a \ a £ G} be independent, 
standard complex Gaussian random variables. Then the eigenvalues {A x | x £ C} of M 
given by f)2. 1 1) are independent, standard complex Gaussian random variables. 

The random matrix ensemble in Proposition 13.11 is the G-circulant analogue of the com- 
plex Ginibre ensemble X, which consists of a square matrix with independent, standard 
complex Gaussian entries. 

Corollary 3.2. Suppose that for each n, \ a S G^} are independent, standard 

complex Gaussian random variables. Then E//( n ) = 7c for each n, and pS- n ^ — > 7c weakly in 
probability. Furthermore, if |G^ n ^| = Q(n £ ) for some e > 0, then — > 7c weakly almost 
surely. 

Proof. For each, say, Lipschitz / : C — > R, 

(E/x)(/) : =E(/i(/)) = ^E/W, 

where the (re) superscripts are omitted for simplicity. By Proposition 13.11 each A x is dis- 
tributed according to 7c, and so (E / u)(/) = 7c(/)- Thus E/i = 7c . 

By the concentration properties of Gaussian measure (see [E]), since the \ x are dis- 
tributed as independent standard complex Gaussian random variables, if / is 1-Lipschitz, 
then 

P[K/)-7t(/)|>t] <2e-l G l* 2 
for each t > 0. If |G^ n ^| = tt(n £ ), then the Borel-Cantelli lemma implies that ^ n \f) — > 
7c(/) almost surely. Applying this to a countable dense family of /, it follows that /i^ — > 7c 
weakly almost surely. 

In the general case, since |G^ n ^| — > 00, each subsequence of /i( n ) has a subsequence 
for which, say, > j, so that by the above argument converges to 7c almost 

surely as j — > 00. It follows that converges to 7c in probability. □ 

The next proposition deals with the G-circulant analogue of the Gaussian Unitary En- 
semble (GUE), which, up to a choice of normalization, is distributed as 2~ 1 ^ 2 (X + X*), 
where X is the complex Ginibre ensemble mentioned above. Equivalently, the diagonal 
entries of the GUE are standard real Gaussian random variables, the off-diagonal entries 
are standard complex Gaussian random variables, and the entries are independent except 
for the constraint that the matrix is Hermitian. It is worth noting explicitly that while each 
entry of the GUE has (complex) variance 1, the variance of a diagonal entry and the real 
part of an off-diagonal entry differ by a factor of 2. (Again, the special case for classical 
circulant matrices was observed earlier in [15].) 

Proposition 3.3. Let G be a finite abelian group and let {Y a \ a £ G} be random variables 
which are independent except for the constraint Y a -i = Y a , and such that 

v Jtk ifa 2 = l, 

Then the eigenvalues {A x | x G ^} of M given by (12. lj) are independent, standard real 
Gaussian random variables. 
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Proof. Let {Z a \ a G G} be independent, standard complex Gaussian random variables. 
Then {Y a \ a G G} are distributed as \2~ x l 2 {Z a + Z a -i) \ a G G}. Thus the eigenvalues 
X x of M in the present proposition are jointly distributed as \/2 times the real parts of the 
eigenvalues of the random matrix defined in Proposition (3TTJ and are thus independent real 
standard normal random variables. □ 

Observe that in the "G-circulant GUE" of Proposition 13.31 every element a G G with 
a = a -1 corresponds to a "diagonal" of M in which the entries are constrained to be real. 

The following corollary follows from Proposition 13.31 in the same way that Corollary 13.21 
follows from Proposition 13.11 

Corollary 3.4. Suppose that for each n, {Ya \ a G C™)} are real and complex Gaussian 
random variables as described in Proposition \3 .3\. Then E//(™) = 7k for each n, and fjy 1 ' —> 
7k weakly in probability. Furthermore, if | | = J7(n e ) for some e > 0, then — > 7k 
weakly almost surely. 

The real Ginibre ensemble X consists of a square matrix with independent, real standard 
Gaussian random variables. The Gaussian Orthogonal Ensemble (GOE) is distributed as 
2~ 1 / 2 (X + X t ). Equivalently, the diagonal entries of the GOE are distributed as N(0, 2) and 
the off-diagonal entries are distributed as N(0, 1). In general the analogues of Propositions 
13.11 and 13.31 for matrices with real entries are less elegant. In the nonsymmetric case the 
eigenvalues have a Gaussian joint distribution in a \G\ -dimensional real subspace of C' G ', 
and in the symmetric case the |G| eigenvalues are not independent in general. We will 
not state such results in general, but will note for future reference that in the "G-circulant 
GOE", every element a G G with a = a -1 corresponds to a diagonal of M in which the 
variance of the entries is 2 instead of 1. (See Theorem 14. 31 below and the discussion following 
it.) 

On the other hand, the analogous results are simple in the case in which the characters 
X G G are all real-valued, so that the Fourier transform defines an isometry (up to scaling) 
between the real I 2 spaces on G and G. By Lemma 12.21 this is the case precisely when 
every a G G satisfies a 2 = 1, or in other words, when G = (^2)"" for some n. In this case 
a G-circulant matrix is automatically symmetric, so that there is no difference (except for 
scaling) between the "G-circulant real Ginibre ensemble" and the "G-circulant GOE" . The 
following results are proved in the same way as Proposition 13.11 and Corollary 13.21 

Proposition 3.5. Let G = (2^)™ and let {Y a \ a G G} be independent, standard real 
Gaussian random variables. Then the eigenvalues |A X | x £ ^} of M given by (|2.1I) are 
independent, standard real Gaussian random variables. 

Corollary 3.6. Suppose that for each n, G^ = (2^)™ and {Ya^ \ a G G^} are indepen- 
dent, standard real Gaussian random variables. Then E/i( n ) = 7^ for each n, and — > 7k 
weakly almost surely. 

4. General matrix entries 

Our main results are stated under a Lindeberg-type condition on the random variables 
Ya^ used to generate the random matrices: 

w Ve > : E E (l^ (n) l 2l i^i> £v ^) = °- 

1 1 aeGW 
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The usual remarks apply about the sufficiency of identical distribution or a Lyapunov- 
type condition: (|4,ip holds in the settings of Theorems 14.11 and 14.21 if all the Y^ are 
identically distributed, or have uniformly bounded (2 + 5) moments; it holds in the settings 
of Theorems 14. 31 and 14. 41 if all the random variables with a given variance assumption satisfy 
such assumptions. 

We now state our main results, deferring the proofs until the end of the section. 

Theorem 4.1. Let a G [0, 1]. Suppose that for each n, {Ya | a G G^} are independent; 
that 

EY™ = 0, E I yj n) | 2 = 1 , and E (Yj n) f = a 

for every a G G^ ; and that (|4.1|) holds. Suppose further that lin^^oo = p exists. Then 
// n ) converges, in mean and in probability, to (1 — p)7c + Via- 

One of the main special cases of interest in Theorem 14.11 is when a = 1, that is, when 
the matrix entries are all real. In that case, the limiting spectral distribution of is 
complex Gaussian if the number of a with a 2 = 1 is negligible for large n. On the other 
hand, if the fraction of such a is asymptotically constant then, due to the presence of many 
real-valued characters x, the limiting spectral distribution will be a mixture of 7c and 7r. 

The other main special case of interest is when a = 0, so that the matrix entries have un- 
corrected real and imaginary parts. In that case, which generalizes the setting of Corollary 
13.21 one can remove the assumption that p^ approaches a limit. 

Theorem 4.2. Suppose that for each n, {Y a \ a G G^} are independent; that 

EYj n > =0, E|Yj n )| 2 = 1, and E(y] n )) 2 = 

for every a G G^ ; and that (|4,ip holds. Then converges, in mean and in probability, 
to 7c- 

The special case of Theorem 14.21 for classical circulant matrices (that is, when the G^ 
are cyclic groups) was proved by the author in |15j . 

Theorem 4.3. Let a G [0, 1], /3 > 0. Suppose that for each n, {Y a (n) | a G are mean 

and independent except for the constraint Y^+ = Ya ! that 

1 if a = b~ 1 ^a- 1 , 
a if a = b ^ a -1 , 
/3 if a = b = a -1 , 
otherwise, 

for a, b G G^ ; and that (|4.ip holds. Assume further that lim n _ > . 0O = P exists. Then 
converges, in mean and in probability, to 

(1 - p)DSf(0, 1 +p(p - a - 1)) + pN(0, 1 + a + p(p - a - 1)). 

ifp < 1 and to N(0,/3) if p = 1. 
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Observe that by Lagrange's theorem on orders of subgroups, l/p\ is an integer, which 
implies that if p < 1 then in fact p < 1/2, and therefore the stated variances of the normal 
distributions named above are indeed positive. 

The most obvious (though not necessarily, as we shall see, the most natural) special case 

of interest in Theorem 14,31 is when the Y^ are real and i.i.d. (except for the symmetry 
constraint), so that a = (3 = 1. In that case the limiting spectral distribution is the mixture 
distribution 

(4.2) (l-p)>T(0,l-p)+pN(0,2-p). 

Two other special cases are suggested by considering the analogy with the GOE and 
GUE. The G-circulant analogue of the GOE, as discussed in the previous section, would 
have real entries such that a = 1 and f3 = 2, and thus the limiting spectral distribution 

(4.3) (l-p)X(0,l)+pN(0,2). 

The slightly simpler nature of this limiting distribution (note that the parameter p plays only 
one role in (|4.3p . as opposed to two roles in (|4.2p ) reflects that a "GOE-like" normalization 
of entries is more natural than equal variances. However, this phenomenon is only evident 
when < p < 1 . In the classical case of Wigner matrices it is well known that in order for 
the semicircle law to hold, no variance assumption need be made on the diagonal entries of 
the matrix. The situation described above emphasizes that this is the case precisely because 
the number of diagonal entries in a Wigner matrix is negligible. 

Finally, when the second moments are the same as for the "G-circulant GUE" of Propo- 
sition [3]3l then a = and (3 = 1 and, as in Corollary 13.4} the limiting spectral distribution 
is simply the standard real Gaussian distribution, even regardless of the value of p. Thus for 
G-circulant matrices, a constraint to be complex Hermitian appears to be somehow more 
natural than a constraint to be real symmetric. As in Theorem 14.21 the assumption that P2 
approaches a limit can even be removed in this situation. 

Theorem 4.4. Suppose that for each n, {Ya \ a £ G^} are mean and independent 

except for the constraint Y^\ = Y a (n) ; that E | yj n) | 2 = 1 for every a e ; that E ( Y a (n) ) 2 = 
if a ytz a -1 ; and that (|4.ip holds. Then /i( n ) converges, in mean and in probability, to tr. 

The special case of Theorem 14.41 for classical circulant matrices (with more restrictive 
assumptions on the distributions of the matrix entries) was proved by Bose and Mitra in 

m- 

We will not attempt to deal thoroughly with the question of when the convergence in 
probability in the results above can be strengthened to almost sure convergence. However, 
the following result gives some sufficient conditions. Each of the conditions stated auto- 
matically implies the Lindeberg-type condition (|4.ip ; for the first part this follows from 
exponential tail decay which is implied by a Poincare inequality (see |12t Corollary 3.2]), 
and for the other parts it is elementary. 

Theorem 4.5. In the setting of Theorem \4-l\ \4-3\ or \4-4\ suppose in addition that 
|G( n )| = f2(ra e ) for some e > and that one of the following conditions holds: 

(1) There is a constant K > such that for every n and every a G G^ n \ Ya satisfies 
a Poincare inequality with constant K. That is, 

Var/(yi"))</,E|v/(yi"))| 2 
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for every smooth f : R 2 — > R. 

(2) There is a constant K > such that \Y^ | < K a.s. for every n and every a 6 G^ n \ 

(3) For some 5 £ (0, 1], sup n6N max s6G („) E|yJ"^| 2+5 < oo, and X^J^Li |G^ n ^| 5 ^ 2 < oo. 

(4) For some S G (0, 1], sup neN max aeG („) E| 1^ | 2+5 < oo, and p^ — > p > 0. 
Then fj,( n ' converges to the stated limit almost surely. 



We now turn to the proofs of our main results. Unsurprisingly, generalizing the results 
of the last section to non-Gaussian matrix entries is achieved by using an appropriate 
version of the central limit theorem to show that the eigenvalues \ x are approximately 
distributed like uncorrelated Gaussian random variables. Even to prove asymptotic results, 
it is necessary here to apply some quantitative version of the central limit theorem, in order 
to achieve suitably uniform control over the X x . The approach taken here (and earlier in 
|15j ) generalizes and extends the method used by Bose and Mitra in [3], which applied 
a multivariate version of the Berry-Esseen theorem and thus required the matrix entries 
to have uniformly bounded third moments. Here a quantitative, multivariate version of 
Lindeberg's theorem is applied. 

If / : M. d — > R is bounded and Lipschitz with Lipschitz constant \f\ L , its bounded 
Lipschitz norm may be defined by 

||/|| Bi = nu«{||/|| 00> |/| i }. 

The bounded Lipschitz distance between random vectors X and Y in R d is defined by 

d BL (X,Y)= sup |E/(X)-E/(Y)|. 
ll/IU<i 

It is well known (see e.g. [HI section 11.3]) that the class of bounded Lipschitz functions 
is a convergence-determining class. The subclass of compactly supported such functions is 
furthermore separable with respect to the sup norm [H Corollary 11.2.5]. Thus to show 
that a sequence of probability measures on R d converges weakly to v in mean, in 
probability, or almost surely, it suffices to show that for each bounded Lipschitz function /, 
v (/) ~~ ^ v {f) in l ne same sense. 

The following is a special case of [21 Theorem 18.1] (cf. the proof of [2j Corollary 18.2]). 

Proposition 4.6. Suppose that X\, . . . ,Xk are independent mean random vectors in M, d 
such that \ Yj)=i Cov PQ) = h- For e > let 

*( £ ) = rX>(ll^Vj>^)- 

3=1 

Then 

d BL f^XX'^j <^ o inf i ( £ + ^)), 
where Z is a standard Gaussian random vector in R d , and > depends only on d. 
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Proof of Theorem\4A\ Let / : C -)• R with ||/|| BL < 1. Observe that 



(4.4) 



Mf) = rk E E /( A x) = E E / ( 7p E x(«)r a 

X6G X 6G a€G 



where (n) superscripts have been omitted for simplicity. We consider X x as a sum of 
independent random vectors in R 2 = C. The relevant covariances are 



Cov( X (a)Y a ) 
The identities 

(4.5) 



E(Re x(a)Y a ) 2 E(Re X (a)Y a )(lm x (a)Y a ) 

E(Re x(a)Y a ) (Im X (a)F a ) E(Im x(a)^a) 2 



(Re w) (Re z) = i Re \(w + w)z], 
(Im«))(Imz) = | Re [(To — 10)2] , 
(Rew;)(Imz) = |lm[(t 



ity — u; z 



will be useful. 

Setting to = z = x( a )Y a for a fixed x £ C, 



EE(Rex(«)^) 2 = E 



aeG 



aGG 



^ReE(x(a) 2 y a 2 + | X (a)| 2 |y a | 2 



I^|G|+aj;x 2 (a))=^(l + al x ^). 



In the last step we have used that unless x 1S real-valued, x an d X are distinct characters, 
and hence orthogonal in £ 2 (G). In similar fashion, we find that 

Cov(A x ) = p E C °v(x(a)Y a ) = \{h + l x =x[S _°«]). 



aeG 



Observe in particular that if a = 1 and x is real-valued, then X x is almost surely real, with 
variance 1; in that case we treat X x as a random variable in R, as opposed to a random 
vector in R 2 . Proposition 14.61 and (|4.1|) (recalling that |x( a )| = 1 always) now imply that 
there is a sequence 5 n decreasing to such that for each x £ G, 



if x is real-valued, and 



E/(^Wl>(a)Yj "-**(/) 

\Vl G l aeG / 

E/(-^^x(a)y a ] - 7 c(/) 



<0~n 



<<5n 
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(4.6) 



x(a)Y a J -p 2 7a(/) 

x(a)Y a ) -(l-p 2 ) 7c (/) 



< P2$n + (1 -P2)S n = 8 n , 

where as above the subscripts (n) are omitted. Since p^ — > p, it follows that =4> 
(1 - p)j c + Via, and so E^( n ) =^ (1 - p)j c + Pla- 
Next observe that 

(4.7) E(M/)) 2 = ^2 E E /( A xi)/( A x,) = E Ei K A xi> A x 2 ))i 

Xi,X2GG XI1X26G 

where F : C 2 — )• M is defined by -F(u>, z) = f(w)f(z), so that H-FH^ < 2. We now consider 
(A Xl ,A X2 ) as a sum of independent random vectors in R 4 . The upper-left and lower-right 
2x2 blocks of Cov((A Xl , A X2 )) are of course just Cov(A Xl ) and Cov(A X2 ), computed above. 
For the off-diagonal blocks, we use w = xii a )Y a and z = X2{o)Y a in (|4.5j) to obtain for 
example 

"1 



Y j n^Xi(a)Y a )(Re X 2(a)Y a ) = £ 

aeG a6G 



-fa 



ReE [xi(a)x2(a)Y a + xi(a)X2(a) \Y a \ 



E Xi(a)X2(a) + E Xi( a )*2( 



a6G 



a6G 



= ( a ^Xi=X2 + lxi=X2) • 

Similarly, it follows that the off-diagonal blocks of Cov((A xl , A X2 )) are unless xi = X2 or 
Xi = X2- 

Assume for now that xi X2 and xi X2- Applying Proposition 14.61 we now obtain 
that there is a sequence 5' n decreasing to such that whenever < 1, 

|E/(A Xl )/(A X2 )- 7a (/) 2 | <S' n 

if xi and X2 are both real-valued, 

|E/(A Xl )/(A X2 )- 7C (/)7«(/)l <S'n 
if exactly one of xi and X2 is real-valued, and 

|E/(A X1 )/(A X2 )- 7C (/) 2 | <S' n 

if neither xi nor X2 is real-valued. (Note that Proposition 14.61 may be applied in the case 
of nonidentity covariance via a linear change of coordinates. For a < 1, the determinant of 
the covariance is bounded away from zero, whereas for a = 1 the variables are real.) Given 
Xi, note that there are at most 2 characters X2 which are unaccounted for. By (|4.7j) . it now 
follows that 

2 



(4.8) 



E M W(/) 2 -i»(/) 2 <S' n + 
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Finally, 



E 



M W (/) _ „(») (/) = [ Em W (/ )2 _ „(») (/) 2] _ 2z » (/) [ E/4 („) (/) _ ^(n) (/) ] 



< 



E/>)(/) 2 -z»(/) 2 +2 -*»(/) 



so by (JO]) and (|Q]> . 



M (n) (/)")• [(l-p)7C+P7a](/) 



in L 2 , and hence in probability. 



□ 



Proof of Theorem \4-%\ The proof is analogous to that of Theorem l4.14 setting a = 0. In that 
case Cov(A x ) no longer depends on whether x is real-valued, which makes it unnecessary 

to assume that approaches a limit. □ 



Proof of Theorem \4-!3\ We omit (n) superscripts as before. We will assume that p < 1; the 
case p = 1 (which implies that in fact pi = 1 for sufficiently large n) is similar and slightly 
simpler. Let A = {a G G \ a = a -1 }. Since G is abelian, A is a subgroup of G. The 
restriction of a character of G to ^4 is a character on A, which is necessarily real-valued on 
A. It follows that for XI1X2 £ G, 



|G|EA X1 A X2 = ^ X i(a)X2(b)EY a Y b 

a,beG 

= ^Xi{a)[(x2(a) + ax2{a))t a ^ a -i + /3x2(a)l a= 



(4.9) 



aeG 

= Xi(a)X2(a) + a ^ Xl(a)X2(a) + /? ^ Xi(o)X2(a) 

aeG\A aeG\A aeA 

= X] Xi(o)x2(o) + « X Xi( a )*2(a) + (/?-«- 1) ^ Xi(a)X2(a) 

a6G aeG a6A 

= l G l (Ai=X2 + «lxi=X2) + 14 (P ~ a ~ ^xik^U 
= \ G \ (lxi=X2 + Ql xi=» +P2(/3-a - l)l X iU=X2u)- 
In particular, for x £ G, 

Var(A x ) = 1 + al x=x +p 2 (/3 - a - 1). 

Denoting 

„(«) = (1 - ^ n) )N(0, 1 + j4 n) (/3 - a - 1)) + ^ n) X(0, 1 + a + (/3 - a - 1)) , 
it follows as in the proof of Theorem gj] that d B L(^H {n \ ^ (n) ) -> 0, and thus that 
E/iW (1 - p)N(0, 1 + p(p - a - 1)) + pN"(0, l + a + p(P-a-l)). 

In this situation just the 1-dimensional case of Proposition 14.61 is necessary. Observe also 
that the variances Var(A x ) are uniformly bounded away from (cf. the comments following 
the statement of the theorem.) This is necessary so that Proposition 14.61 may be applied 
for nonidentity covariance, via a linear change of coordinates, and still yield error bounds 
5 n which are uniform in / with < 1. 
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By ((HI), if xi / X2 and xi X2> then 

Cov((A Xl ,A X2 )) = (l+p 2 (/3-a-l))l 2 + a 
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-"-Xi=Xi 

1 





X2=X2. 



+ p 2 (P - a - l)l xl]A=X2lj 



1 

1 



We consider separately the cases p = and p > 0. If p = 0, then when xi 7^ X2 and 
Xi X2, we have 



Cov((A Xl ,A X2 )) 



1 + al Xl=xT 

1 + al X2=x5 



From here the argument is completed as in the proof of Theorem 14.11 
Suppose now that p > 0. Given ) 

X2 € G with xx\a = X2\A- Therefore, 



Suppose now that p > 0. Given xi £ G, by Lemma 12.31 there are exactly — values of 



Cov((A Xl ,A X2 )) = (l+p 2 (/3-a-l))/ 2 + 



lx2=X2. 



for all but a negligible fraction of pairs Xl)X2 £ G. The argument is again completed as in 
the proof of Theorem 14.11 □ 



Proof of Theorem \4-4\ The proof is analogous to that of Theorem 14.31 setting a = and 



(3 = 1. In that case Cov(A Xl ,A X2 } 
approaches a limit. 



1 



Xl=X2' 



so it is unnecessary to assume that p\ 



(n) 



□ 



Proof of Theorem \4-5\ (1) The Poincare inequality assumption and independence imply 
an exponential concentration property for the family of eigenvalues {A x | x G G^}. 
In particular, combining Corollaries 5.7 and 3.2 of [12], it follows that for each L- 
Lipschitz F : £ 2 (G^) -> R, 

F(yW) - EF(YW) > t] < 2e- ct/VFL 
for each t > 0, where c > is some absolute constant and y( n ) is shorthand for 



aeG(") - 



Now for a 1-Lipschitz / : C -)• R and A; £ N 
A- 



i=i 



< 



UU — 2,' 



by the Cauchy-Schwarz inequality. Combining this with Lemma I2.1l f2"j) it follows 

— 1 /2 

that M (n) (/) is |G (n) | 1 -Lipschitz as a function of F" , and so 



//»>(/) -E/>)(/)) 



> t 



< 2e 



-ctJ \/K 



Combined with the already known convergence in mean and the Borel-Cantelli 
lemma, this implies almost sure convergence of jJ-(f). 
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(2) The proof is similar to the previous part, using instead Talagrand's convex-distance 
concentration inequality for independent bounded random variables \17\ Theorem 
4.1.1] (see e.g. |141 Corollary 4] for an explicit statement of a version that applies 
directly to complex random variables), cf. the proof of [15\ Theorem 2]). 

(3) The stated Lyapunov-type assumption yields upper bounds on all the 5 n quantities 

in the proofs above of order |G( n )| for < 5 < 1 (cf. [2J Corollary 18.3]). Thus 

the assumption that Yln°=i \ | allows the Borel-Cantelli lemma to be applied 
again. 

(4) The assumption that p > implies that |G^ n ^| actually grows exponentially: since 

p^ is always the reciprocal of an integer (by Lagrange's theorem about the orders 

of subgroups of finite groups), p^ — > p > implies that i s eventually constant. 
By the classification of finite abelian groups, 

where m > 0, kj > 1 for each j, and each nonidentity element of H has odd 
order. (For simplicity of notation, we are again suppressing the dependence of all 
these on n.) In this notation, the number of a £ G such that a = a -1 is 2 m , so 
that |£r| = 2 m /p2- The hypothesis that |G^ n ^| is strictly increasing thus implies 
that m is eventually strictly increasing, and hence |G^ n ^| is eventually exponentially 
increasing. Therefore the previous part of the theorem applies. □ 



References 

T. M. Apostol. Introduction to Analytic Number Theory. Undergraduate Texts in Mathematics. 
Springer- Verlag, New York, 1976. 

R. N. Bhattacharya and R. Ranga Rao. Normal Approximation and Asymptotic Expansions. Robert E. 
Krieger Publishing Co. Inc., Melbourne, FL, 1986. Reprint ol the 1976 original. 

A. Bose and J. Mitra. Limiting spectral distribution ol a special circulant. Statist. Probab. Lett., 
60(1):111-120, 2002. 

A. Bose and A. Sen. Another look at the moment method lor large dimensional random matrices. 
Electron. J. Probab., 13:no. 21, 588-628, 2008. 

A. Bose, R. Subhra Hazra, and K. Saha. Spectral norm of circulant type matrices with heavy tailed 
entries. Electron. Commun. Probab., 15:299-313, 2010. 

W. Bryc and S. Sethuraman. A remark on the maximum eigenvalue for circulant matrices. In High 
Dimensional Probability V: The Luminy Volume, volume 5 of IMS Collections, pages 179-184. Institute 
of Mathematical Statistics, Beachwood, OH, 2009. 

P. Diaconis. Group Representations in Probability and Statistics. Institute of Mathematical Statistics 

Lecture Notes — Monograph Series, 11. Institute of Mathematical Statistics, Hayward, CA, 1988. 

P. Diaconis. Patterned matrices. In Matrix Theory and Applications (Phoenix, AZ, 1989), volume 40 of 

Proc. Sympos. Appl. Math., pages 37-58. Amer. Math. Soc, Providence, RI, 1990. 

R. M. Dudley. Real Analysis and Probability, volume 74 of Cambridge Studies in Advanced Mathematics. 

Cambridge University Press, Cambridge, 2002. Revised reprint of the 1989 original. 

S. R. Jain and S. C. L. Srivastava. Random cyclic matrices. Phys. Rev. E, 78:036213, 2008. 

J. Kahane. Some Random Series of Functions, volume 5 of Cambridge Studies in Advanced Mathematics. 

Cambridge University Press, Cambridge, second edition, 1985. 

M. Ledoux. The Concentration of Measure Phenomenon, volume 89 of Mathematical Surveys and Mono- 
graphs. American Mathematical Society, Providence, RI, 2001. 

M. B. Marcus and G. Pisier. Random Fourier Series with Applications to Harmonic Analysis, volume 
101 of Annals of Mathematics Studies. Princeton University Press, Princeton, N.J., 1981. 



THE SPECTRA OF RANDOM ABELIAN G-CIRCULANT MATRICES 



15 



[14] M. Meckes. Concentration of norms and eigenvalues of random matrices. J. Fund. Anal., 211(2):508- 
524, 2004. 

[15] M. Meckes. Some results on random circulant matrices. In High Dimensional Probability V: The Luminy 
Volume, volume 5 of IMS Collections, pages 213-223. Institute of Mathematical Statistics, Beachwood, 
OH, 2009. 

[16] J. Serre. Linear Representations of Finite Groups, volume 42 of Graduate Texts in Mathematics. 

Springer- Verlag, New York, 1977. Translated from the second French edition by Leonard L. Scott. 
[17] M. Talagrand. Concentration of measure and isoperimetric inequalities in product spaces. Inst. Hautes 

Etudes Sci. Publ. Math., (81):73-205, 1995. 
[18] W. Yin, S. Morgan, J. Yang, and Y. Zhang. Practical compressive sensing with Toeplitz and circulant 

matrices. Rice University CAAM Technical Report TR10-01, 2010. 

Department of Mathematics, Case Western Reserve University, 10900 Euclid Ave., Cleve- 
land, Ohio 44106, U.S.A. 

E-mail address: mark.meckes@case.edu 
URL: www. case . edu/artsci/math/mwmeckes/ 



